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چکیده در اين مقاله» به شرح روشی برای بهبود میدان تنش به دست آمده از روش ایزوژئومتریک و برآورد حطای موجود در آن پرداحته 
می‌شود. این ب رآورد کننده نعطا ( 65111۳16010۲ 12۲۲0۲) در زمره دسته روش‌های برآورد عطای مبتنی بر بازیافت تنش قرار می‌گیرد. در این روش, با 
استفاده از تقطه‌های فراهمگرا؛ برای تابع مفدارهای هریک از مولفه‌های میدان تنش در هر ناحیه. یک سطح فرضی ساحته می‌شود. برای تعریف این 
سطح همان تابع‌های شکل نریزی (01610۳75] 5/10۳6 "حفلل]۷) به کار می‌روند که با کمک آن‌ها در روش ایزوژئومتریک تابع جابجایی تقریب 
زده می شود. از مفایسه ثرم عطای کارمایه دقیق و رم حطای کارمایه تقریبی, مشاهده می‌شود که برآورد کننده حعطای پیشنهادی ا زکارایی مناسبی 


برای ب رآورد نعطای موجود در تحلیل مساله‌ها به روش ایزوژئومتریک برنحوردار است. 


واژه‌های کلیدی برآورد خطا بازیافت تنش, تحلیل ایزوژئومتریک 
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۳ اول مقاله در تاریخ ۷۱ و سخه‌ی نهایی آن در تاریخ ۳ به دفتر نشربه رسیده استتت: 
)۱ نویسنده‌ی مسوول. دانساز؛ دانشکده عمران. فاسش‌گاه صنعتی شاهرود 

(۲) دانشجوی دکترای سازه دانشکده عمران دانشگاه صنعتی شاهرود 

(۳) استادیا دانشکده عمران دانشگاه صنعتی شاهرود 


۱/۸ 


مقد مه 

با پیشرفت شتابان دانش و فن آوری» روش‌های عددی 
گسترش و گوناگونی چشمگیری یافته‌اند. از جدیدترین 
این روش‌هاء می‌توان روش ایزوژئومتریک را نام برد. در 
روش تحلیل ایزوژئومتریک. ویژگی‌های بی‌مانند و 
مناسبی نهفته است که شاید در آینده‌ای نه چندان دور 
بتواند جایگزین روش‌های عددی معمول نظیر اجزای 
محدود و روش‌های بدون شبکه (60000ظ 1۲66-ظوع۷۲) 
شود. وجود برخحی کاستی‌ها در روش اجزای محدود. 
مانند وجود خطا در تعریف مرزهای مسأله‌های با هندسه 
پیچیده و یا مسأله‌های با تغییرات شدید در بارگذاری و 
و کی هاش ما وت مت لا تانق تهج 
در برخی مسأله‌ها؛ همانند مسأله‌هایی که در چارچوب 
لاگرانژی حل می‌شوند و يا مسأله‌های بهینه‌سازی شکل 
سازه. از دسته علت‌هایی هستند که به نوآوری این روش 
انجامید. این روش را برای نخستین بار هیوز و 
همکارانش در سال ۲۰۰۵ معرفی کردند [1]. روش 
ایزوژئومتریک از برخی مفاهیم همانند روش اجزای 
محدود و روش‌های بدون شبکه برخوردار است که 
می‌تواند برای حل مسأله‌های معادلات دیفرانسیل با 
مشتقات جزئی حاکم در زمینه‌های گوناگون علوم و 
مهندسی, از جمله نظریه کشسانی» به کار رود. در این 
کب کب 0 ی 
(معنعع(1 ۸060 ۲عاتا0۳))) بهره گرفته شده است. 
دلیل اینکه روش‌های طراحی به کمک رایانه تاکنون در 
روش اجزای محدود به صورت گسترده وارد نشده است 
به اختلاف زمانی» در پیدایش این دو نسبت به یکدیگر 
بستگی دارد. آغاز پیدایش روش اجرای محدود در سال 
ادف بوده اش هر ای که روش قاع ط نف 
به کمک رایانه بعدها در حدود سال‌های ۱۹۷۰ تا ۱۹۸۰ 
شکل گرفته‌اند [1] 

انديشه بنیان گذاری روش ایزوژئومتریک بر پایه 
بی - اسپلاین‌های نسبی غیر یکنواخت «08-1/011018) 
۱-0 210821ع) (9ظل1ل) می‌باشد. در این روش 


نشریه مهندسی عمران دانشگاه فردوسی مشهد 


برآورد نحطا وه مان کر یوکس یهن و 


با کاربرد ویژگی‌های تابع‌های پایه اسپلاین و نربز در 
تعریف دقیق خم‌ها و سطح‌هاء از آن‌ها برای درونیابی و 
تقریب‌سازی هم استفاده می‌شود. باید یادآوری شود با 
آنکه شرح دادن روش ایزوژئومتریک مورد نظر نمی‌باشد 
ولی برای آشنایی خواننده با مفهوم‌هایی که در برآورد 
خطای این روش با آن‌ها روبرو خواهيم شد. در 
بخش‌های بعدی به بیان بعضی از مفهوم‌های نخستین 
پداته له آستتا ام | با هی تن مشرم وان ی 
مرجم‌های [1-7] مراجعه کرد. 

خطاء بخش جدانشدنی تحلیل‌های عددی به شمار 
می‌رود و همواره پژوهشگران را در قابلیت اعتماد نتيجه‌ها 
نگران می کند. در حالت کلی» روش‌های برآورد خحطا در 
دو دسته روش‌های بازيافت تسش و روش‌های 
باقیمانده‌ای قرار می‌گیرند [8]. در پژوهش حاضر یک 
روش هندسی وابسته به بازیافت تنش. و همسو باروش 
ایزوژئومتریک» برای برآورد حطای نتیجه‌ها پيشنهاد شده 
است که همه مرحله‌های آن در بخش‌های بعدی شرح 


داده خحواهد شبال. 


بی - اسپلاین و نربز 

این بخش در حد نیا به معرفی کوتاهی از خم‌هاو 
سطح‌های بی- اسپلاین و نربز می‌پردازد. برای آشنایی 
بیشت مراجعه به مرجع‌های [9,10] پيشنهاد می‌شود. 

بر ما ار 
بی -اسپلاین ها در یک فضای پارامتری (ناحیه با ط6ا۳۵) 
تغریف می‌شوند, این ناحیه‌ها دامن الگوسازی شده وا بنه 
جندین زیر دامنه تقسیم می‌کنند. یک بردار گرهی « :۹۹ 
۲ در فضای پارامتری یکبعدی از یک دنباله 
مختصات به صورت زیر تشکیل می‌شود[10]. 


)۱( 0 رم وصول ک( 2 ببه : 9 
که در آن بیق امین گره, 0 مرنبه چند جمله‌ای و 9 شمار 


تابع‌های شعل تشکیل‌دهنده یی اسپلاین به شمار 
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می‌رود. چندین گونه از بردارهای گره‌ای وجود دارد ولی 
در این بحث تنها از گونه ویژه‌ای از بردارهای گره‌ای» به 
نام بردارهای گره‌ای نامتناوب ( ءمصع1 عذ8منتهم‌مما 
۲ (یا باز (0060)) بهره گرفته می‌شود. این گونه 
بردارها به شکل زیر نشان داده:من شوند: 


)۲( مرگ و - > 


2+1 241 


درحه «ر که با (ج) و تسیا داده می‌شود به صورت 


زیر تعریف می‌شود [10: 


ی 11 1 
۳ 1+1 1 » ۷ 
1 روت 9:9۱ ۲ 012 
بمب وت 
ات ی و ایا ان اه 
كِِ اه 1 و 9 


درجه 9 به صورت زیر پدید می‌آید[10]: 
() «اع>و > ه (۷:,)6( 2 -(0)6 


(ع)) یک خم چند جمله‌ای قطعه ای 
(وسه امتجمصرامم ۳۱6۵۵۷/[56) است که در آن ۳ 
نقطه‌های هادی و ((ع),,(] تابع‌های پایه‌ای بی- 
اسپلاین هستند که روی بردار گره‌ای نامتناوبی به صورت 
رابطه (۲) با فرض 2-0 و 9-1 تعریف می‌شوند. 
اگر 0 درجه تابع‌های پایه. 1+ شمار نقطه‌های هادی و 
1+ شمار گره‌ها باشند آنگاه می‌توان رابطه 
041 +20 0 را برای آن‌ها نوشت. 

با روش مشابهی سطح‌های بی- اسپلاین با رابطه 


(۵) تعریف می‌شوند|[10]: 
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1 
)0( زط() پل( (ع) مر - ٩),1(‏ 
0-ز 1-0 
که در آن: 
0( : سک وی > 
941 8+1 
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به گونه‌ای که بردار گره ای 2 دارای ۲+1 گره 


و دارای 5+1 گره می‌باشد . 


(5) ( 
0 9 >و > ه 2( 
۷() :2 


که در آن ۴۱ نقطه‌های هادی ۷,۱) وزن‌ها 
و ((ع) ,,1) تابع‌های پایه‌ای بی- اسپلاین از درجه 8 
هستند. که بر روی بردار گره‌ای به صورت رابطه(۲) 
تعریف شده‌اند. 

و در پایان. یک سطح نربز که در جهت و از 
درجه 0 و در جهت 1 از درجه 4 باشد با رابطه زیر 


زر (0),(),( 2 
( اعس>0 ______- ٩‏ 
,۷ (10) بر (ج) ,رمرم 


در عبارت بالا |ر:2) شبکه نقطه‌های هادی 
می‌باشد که در دو جهت تعریف شده است؛ همچنین 
[ر,۳] وزن‌هاو ((ع) ,لا و (()پرلا] تابع‌های 
پایه‌ای ی ی اسپلاین هستند که بر روی بردارهای گره‌ای 
با رابطه )۱( تعریف شده‌اند . در رابطه )۸ ا کت تابع‌های 


یایه‌ای نسبی قطعه‌ای مانند رابطه )٩(‏ تعریف شود: 
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۵۱ ,۵( ووزا (ج) ,بل - (1:) 9۱ 


۳ 


رابطه (۲۰) برقرار است: 


[ 


( یرت 3 < :)۹ 


حاصل از آن با تابع‌های پایه درجه دو در جهت : و۷ 
مشاهده می‌شود. 


الف - شبکه نقطه‌های هادی 


۱ 
۱ 


رم 


«ِ 


ب- سطح نربز 


شکل ۱ شبکه نقطه‌های هادی و سطح نربز مربوط به آن با تابع‌های 


پایه درجه دو [10] 


روش‌های برآورد خطا در احزای محدود 
برای آشنابی با برحی از روش‌های برآورد خحطای 
تحلیل‌های عددی در این بخش به معرفی شماری از این 
روش‌ها در برآورد خطای تحلیل اجزای محدود پرداخته 


0 
ت 


ی 2 
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برآورد نحطا وه مان کر یوکس یهن و 


از نخستین مقاله‌هایی که در آن یک روش عمومی 
از برآورد خطا بیان شد. مقاله‌هایی است که ریچاردسون 
در سال ۱۹۱۰ نوشته است. روش ریچاردسون در برآورد 
خطا از این حقیقت پیروی می‌کند که خطا در هر طرحی 
از تفاوت محدود. در حالت‌های معمول به اندازه 
شبکه‌های به کار رفته بستگی دارد [8]. کار اصلی در 
برآورد خطارا در سال ۱۹۷۸ بابوشکا (13005 .1) و 
رینبولت ۳61000100 .) آغاز کردند. روش آن‌ها که 
روش باقیمانده‌ای نام گرفت بر این اساس بود که با 
قراردادن حل حاصل از روش اجزای محدود در معادله 
دیفرانسیل اولیه و محاسبه مانده آن و نیز محاسبه پرش 
میدان گرادیان در لبه اجزا؛ برآوردی از خطای حل 
اجزای محدود به دست می‌آید [8]. دستة دیگر از 
روش‌های برآورد خطاء روش‌های بازیافت تنش می باشد. 
به طور کلی بازیافت تنش» روشی برای بالا بردن دقت. و 
پیوسته و هموار کردن میدان تنش با گرادیان میدان 
حاصل از حل اجزای محدود است. در این روش با بهره 
جویی از حل روش اجزای محدود. یک پاسخ نزدیک به 
حل دقیق یا تحلیلی در هر گره یا در نقطه‌ها دلخواه روی 
دامنه محاسبه شده که دقت بالاتری نسبت به حل اولية 
اجزای محدود داشته است. این میدان تنش بهبود یافته 
با این رابطه تعریف می‌شود: 
(6۱ 6( - و 

ان وهای که نی ۱۲ 
تابع‌های شکل مورد استفاده در اجزاست. با کاربرد این 
میدان بهبود یافته. خطای بازیافت به این صورت محاسبه 
3 
۲( 6 - 6 ره 

کته:فر هار۱ 9 لش ای او آنختر ای 
محدود می‌باشد [11]. روش‌های گوناگونی برای بازیافت 
تنش از حل اجزای محدود وجود دارد که از جمله 
ابتدایی‌ترین آن‌ها می‌توان به روش میانگین‌گیری که 
هینتن و کمپیل (020000611 220 1110100) در سال ۱۹۷۶ 
[12] و روش تصویر ما (۱۷۲6۵0 ۲:0[6600۳ ما ) که 
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اودن و براجلی (۳۵۵6۳11 200 0060) در سال ۱۹۷۱ 
۱۵1 هه کیان بر دا رازه کر دب او سسال ۳ وفره 
بازیافت تنش را بر پایه نقطه‌های فراهم گرا ٩۳۴‏ 
۳۳ پچ و زو 
نوآوری کردند و گام بسیار بلندی در بازیافت تنش 
۱ ام وی ری هت ار کی ار 
بهترین و اثر گذارترین روش‌هابرای ؛ برآورد خطا در 
مسأله‌های مهندسی به کار می‌رود. اساس ایین روش بر 
پایه کاربرد نقطه‌هایی به نام نقطه‌های فراهم گرا است که 
دز آنها تایه تست امله او تلیل تفر سین تستیت باه 
ساير نقطه‌ها از دقت بیشتری برخوردار می باشد. درایین 
روش با برازش یک میدان به صورت چند جمله‌ای با 
ضریب‌های نامعین بر روی گرادیان حاصل از روش 
اجزای محدود روی گروه اجزای پیوندخورده به هر 
گره(۳۵/0) میدان گرادیان بهبود یافته تعیین می‌شود. 
پس از آن نیز فن‌های بسیار دیگری را پژوهشگران 
نوآوری کردند که در هر یک به نوعی به بهبود روش 
5 پرداخته می شود که به عنوان نمونه می توان به 
روش‌های ۲۳۳ [15] و ٩۳3۳۴7‏ [16] اشاره کرد. 


شرح روش بازیافت تنش 
به طور کلی در این روش. میدان تنش بهبود یافته برای 
هر مولفة تنش در هر ناحیه به صورت یک سطح فرضی 
در نظر گرفته می‌شود. اين سطح فرضی با کاربرد تابع‌ها 
شکل نربزی که برای گمانه زدن تابع مجهول (جابجایی) 
استفاده شده‌انده به دست می‌آید. یک سطح نربز زمانی به 
دست خواهد آمد که مختصات نقطه‌های هادی آن 
مشخص شده باشد. با تعریف مختصات :و۲ هر نقطه 
هادی با کمک کاربن تنها مولفه مجهول برای تعیین سطح 
بهبود یافتة تتش. مولفه 7 نقطه هادی می‌باشد. روش 
محاسبه مختصات 2 نقطه‌های هادی به گونه‌ای است که 
سطح تنش جدید به دست آمده نسبت به سطح تنش 
قبلی که از تحلیل ایزوژئومتریک حاصل شده است. به 
یک سطح تنش بهبود يافته تبدیل می‌شود. 
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اساس محاسبه مختصات نقطه‌های هادی و د 
وی زوس رت ی برگرفته از 
قی کین نلطه‌هانی آشت کتاور آن‌ها شش خاضل از تلا 
تقریبی نسبت به دیگر نقطه‌ها از دقت بیشتری برخوردار 
می‌باشد. در این نقطه‌ها مرتبه همگرایی شیب یک تابع. 
یک مرتبه از مقداری که از تقریب تابع شکل وابسته به 
حل تقریبی انتظار می‌رود بالاتر است. به همین دلیل این 
هناهام کر هی کاس تفن که سوه ار 
آن‌ها را پارلو مطرح کرد [11]. 

در تحلیل ایزوژئومتریک دو بعدی که دامنه مسأله 
نقطه‌های گرهی به اجزای چهار ضلعی 
تقسیم می‌شود این نقطه‌های فراهمگرا بر نقطه‌های 
گوسی جزء قرارمی‌گیرند. مختصات نقطه‌ها هادی سطح 
بهینه تنش با کمینه کردن فاصله بین این سطح فرضی و 
سطح تنش به دست آمده از حل ایزوژئومتریک در 
نقطه‌های گوس اجزای هر ناحیه با کاربرد روش کمترین 
مجموع مربعات محاسبه می‌شود. 

اگر سطح بازیافتی (بهینه) هریک از ملفه‌های 
راز تن ای هشونا رسیم سا ها 
شکل نربز می‌توان این سطح را در داخل هر ناحیه به 
صورت زير بیان کرد: 

2 ۳ ۷, ۷( ۱۳ 


احز اح1 


که فقو ان شمان تقظ‌ها هادع هن شعفت ۳ و 6 


با توجه به ارایش 


شمار نقطه‌های هادی در جهت ‏ هر ناحیه. ‏ تابع‌های 
شکل نربز و ۳ مختصات نقطه‌های هادی وابسته به صفحه 
تمیق پا شاه | کر 16 و ۳ به ترتیب. بردار تابع‌های شکل 
نربز و بردار مختصات نقطه‌های هادی. به صورت 
رابطه‌های (۱۶) و (۱8) تعریف شود. رابطه (1۳) را 


می‌توان به صورت (۱۱) نوشت. 
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تشریه مهندسی عمرال دانشگاه فردوسی مشهد 


۳۲ 


همان گونه که مشاهده می‌شود. تنها عامل مجهول 
برای تعیین این سطحء مختصات 2 نقطه‌های هادی (بردار 
۳) می‌باشد. برای تعیین این مقدارها باید مجموع 
را | نس نماض را قعاسا: 
ایزوژئومتریک و تنش بهبود یافته را در نقطه‌های گوس 
کمینه کرد. برای این هدف. تابع (۳)۳ به این صورت 
ت ی و2 : 


یک رک 
2 6 2۵/6 (۳)۳ 


از از 
ادن ان تفای اعطاه خمنس | مه از تحلیل 
ایزوژئومتریک و بو با به ترتیب شمار نقطه‌های گوس 


در جهت‌های ‏ و۷ موجود در هر ناحبه می‌باشد. با 


۳)۳( << )8۳۳ 6۳ (۱۸) 
1-1 


که در آن > شمار نقطه‌های گوس موجود در هر 
ناحیه می‌باشد. در پایان با مشتق‌گیری از تابع (108 نسبت 
به مولفه‌های 2 نقطه‌های هادی و برابر صفر قرار دادن آن, 
مختصات نقطه‌های هادی صفحه تنش بهبود یافته به 


(۱4 عم ح و مراد 0 


که در آن: 


1 1 
هت رخجط ی یر 


با داشتن مختصات نقطه‌های هادی. سطح تشون 
و اد( نیز به دست می‌آید. در ادامه. نشان داده 
می‌شود که اين میدان ملفه تتش نسبت به سطح تنش 
حاصل از روش ایزوژئومتریک دقیق‌تر می‌باشد و از 
این‌رو می‌تواند به عنوان یک برآورد کننده توانای خطا 


برای تحلیل ایزوژئومتریک به کار رود. این براورد ۹ 


نشریه مهندسی عمران دانشگاه فردوسی مشهد 


برآورد نحطا وه مان کر یوکس یهن و 


خحطا با در نظر گرفتن اختلاف بین سطح تنش بازیافتی و 
سطح تنش به دست آمده از تحلیل ایزوژئومتریک برای 
هر جزی به صورت تفریبی به معیاری برای تعیین میزان 
خطای موحود در آن جزء دست پیدا کرد. 

در ادامه برای بررسی کار این شتر اوزد سل 
خطاء به مقایسه میانگین خحطای کارمایه تقریبی و میانگین 
خطای کارمابه دقیق برای ی نمونه بنام در نگره 
کشسانی که بیشتر برای بررسی کاراین براوردکننده‌های 


خطا به کار می‌روند» پرداخته شده اه 


میانگین خطای کارمایه 
از معیارهای مختلفی برای تعیین میزان خطابهره 
می‌گیرند. یکی از بنام‌ترین معیارهای بیان خطا معیار 
خطای کارمایه است. بنابه تعریف. میانگین خطای کارمایه 
دقیق تنش برای یک جزء به صورت زیر بیان می‌شود 
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در این رابطه » مقدار دقیق بردار تنش» 6 تنش به 
یگ اهاز تفا مرس 1۳ رشن کشا نی ٩۳‏ 
دامنه جزء می‌باشد. باید دانست. هميشه حل دقیق 
فساله‌های رنگره کشسانی دوردسترتن یتقو نها برای 
برخی از موردهای ویژه. پاسخ نظری وجود دارد. بنابراین 
به جای کاربرد میزان دقیق تنش, میزان بهبود یافته آن 
برای محاسبه میانگین خطای کارمایه به کار می‌رود. 
میانگین خطای کارمایه تقریبی به صورت زیر تعربف 
و 


1 
2 
 ۱(‏ امه وه 3۳06 60| |-|5| ]۰ 
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که در اینجا » تنش بازیافتی و 6 تنش حاصل از 
تحلیل ایزوژئومتریک می‌باشد. در پایان مجموع میانگین 
خحطای کارمایه اجزاء میانگین خطای کارمایه کل دامنه را 
تشکیل می‌دهد. 

در ادامه به ارائه نتیجه‌های تحلیل آمده برای سه 


مسئله نمونه پرداخته می‌شود. 


تیر طره تیموشنکو 
در این بخش نتیجه‌های الگوسازی یک تير کشسان خطی 
همسانگرد در شرایط تننش مستوی باروش 
ایزوژئومتریک و بازیابی تنش‌های آن ارائه می‌شود 
(شکل (۲)). عامل‌های به کار رفته در این تحلیل از این 


قر ارند: 


5-<- ۷ , ۳-1500 , ۳<300 , 2-2 , 10<,ر] 


شکل ۲ تیر طره در شرایط تنش مستوی 


تنش‌های دقیق این مسأله را با رابطه‌های تحلیلی؛ 


۳ 0 ۱ 
۳۲( 2 بر وه : 6-0 : ۱ ۵ 


"۳ 3 
هه بد-1 قسی باه بآ الکوس ری 


تحلیل این تير به روش ایزوژئومتریک از ۱۰۵ نقطه هادی 
و دو ناحبه همانند. استفاده شده است (شعل ۹0 
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ناحیه دوم ناحیه اول 


شکل ۳ نقطه‌های هادی مورد استفاده در الگوسازی تیر طره 


برای نمایش بهتری از کارایی برآورد کننده خحطای 
پیشنهادی با توجه به مرجع [18] از تابع‌های شکل نربز 
مرتبه دو در هر ناحیه به کار برده شده‌اند. بردارهای 


گرهی در جهت‌های 1و وّبه صورت زیر اند. 
۲۳( و [0,0,0.3,0.5,0.7,1,1) < 1 


۱ 0.2,)(.3,)(.4,)(.5,)(.6,)(.7,)(.8,)0.9,1,1), 1 0,0,0 حِ 


با توجه به مرتبه تابع‌های شکل. در تحلیل 
ایزوژئومتریک شمار نقطه‌های مورد نیاز برای تابع اولیه 
گیری عددی به روش گوس از همانند آن در اجزای 
محدود بیشتر است [119. از این رو در این نمونه از نه 
نقطه گوس برای تابع اولیه گیری عددی و نقطه‌های بهینه 
تنش در هر جزء استفاده شده است. 

نمودار تغییرات میانگین خطای کارمایه تقریسی و 
میانگین خطای کارمایه دقیق برای هر جزء در شکل (1) 
آمده است. بنابه شکل (:. همانندی در چگونگی 
تغییرات میانگین خطای تقریبی نسبت به میانگین خحطای 
دقیق نشانه کارایی مناسب محاسبه‌گر خحطای پیشنهادی 
له 

همچنین شاخحص اثر گذاری برای کل دامنه ۰/۸۱ 
محاسبه شده است. شاخص ار گذاری نسبت مجموع 
میانگین خطای تقریبی به میانگین خطای دقیق می‌باشد که 


نشانگر میزان همگرایی پاسخ به سوی پاسخ واقعی است 


تشریه مهندسی عمرال دانشگاه فردوسی مشهد 


۷ برآورد عطا و بهبود میدان تتش به‌دس ت آمده از ... 


و معیاری برای بیان دقفت محاسبه‌ گر خطابه شمار پیدا کند [11]. در شکل (۵) چگونگی توزیع میانگین 
می‌رود. هنگامی یک محاسبه‌گر خطا دارای کارایی خحطای کارمایه دقیق و میانگین خطای کارمایه تقریسبی 
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شکل ۶ نمودار تغییرات میانگین خطای کارمایه تقریبی و میانگین خطای کارمایه دقیق تیر طره تیموشنکو 


1 ۰-۰" 1( 70 18 ظ 40 34 2 10 


الف- میانگین خطای دقیق 


5 (1 ۰ ۰ 3۵ 1 0 98 


ب- میانگین خطای تقریبی 


شکل ۵ توزیع میانگین خطای دقیق و میانگین خطای تقریبی تیر طره تیموشنکو 


صفحه نامحدود سوراخدار مساله دو بعدی دیگری که برای مقایسه مورد بررسی قرار 
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(شکل ()). 
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شکل 1 صفحه نامحدود سوراخدار 
شده است (شکل (0۷). 

سم 
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شکل ۷ ویژگی‌ها و شرایط مرزی دامنه الگوسازی شدهٌ صفحه 


نامحدود سوراخدار 


به ان صفحه تحت تنش کششی 21 ,» وارد 
شده و با فرض شرایط تنش مسطح تحلیل شده است. 
ویژگی‌های مصالح به صورت کشسان خطی همسانگرد با 
ضریب کشسانی ۲-1000 و ضریب پواسون ۷۲۶0.3 در 
نظر گرفقهع شنده اسنت: تسی‌هانن دقیق. این مساله با ترجه 
به حل تحلیلی آن در روابط (۲۶) تا (۲7) آمده‌اند [17]. 


سال بیست و دوم شماره دی ۱۳۹۰ 


۳۵ 


2 4 


۲( 0 وس 2 + (40 ومه + 20و60 )1-3 ری 
۲ 2 2 ۲ 


4 


2 
97 - - 3-0920 - 009 40( - 23-0 ۲۵( 
۲" 2 1 


/ 


۵ 3 ار 3 
)۳۹( 0 2-10 - + (40 وزو + 20 910-)_ت - - 
ِ ۳ ۳ 1 


6 


سوراخدار به روش ایزوژئومتریک از ۱۳ نقطه هادی و دو 
ناحبه همانند استفاده شده است (شکل ری 


۰ ۱ 
تاحیه اول 


تاحیه دوم 


شکل ۸ نقطه‌های هادی به کار رفته در الگوسازی صفحه نامحدود 


سوراخدار 


نمونه اول. در این نمونه نیز تابع‌های شکل نربز 
مرتبه دو در هر ناحیه و نه نقطه گوسی به کار می‌روند. 
بردارهای گرهی در جهت‌های و چّ برای هر ناحیه به 


صورت زیر در نظر گرفته شده اند. 


1[  )0,0,0.3,0.7,1,1! و‎ ۳۷ 
ج‎  )0,0,0.2,0,3,0.4,0۵.5,0.6,0.7,0.8,1,1[ 


نمودار تغییرات میانگین خطای کارمایه تقریسی و 
میانگین خطای کارمایه دقیق برای هر جزء در شکل )٩(‏ 
به نمایش فرافازی ات 

با توجه به عدد ۰/۸۱ برای شاخحص ار گذاری 
میانگین خطای کل دامنه و نیز همسانی تقریسی نمودار 
تغییرات میانگین خطای کارمایه دقیق و تقریبی در شکل 


)٩(‏ می‌توان نتیجه گرفت که روش برآورد خحطای به کار 


تشریه مهندسی عمرال دانشگاه فردوسی مشهد 


۳۹ برآورد نحطا و بهبود میدان تنش به‌دست آمده از ... 
گرفته شده برای صفحة نامحدود سوراخدار نیز به نسبت میانگین خحطای کارمایه و همچنین مقایسه آن بین میانگین 


از کارایی خوبی برخوردار می‌باشد. خطای تقریبی و دقیق, به رسم توزیع میانگین خطای 
در شکل (۱۰) برای دریافت بهتری از تغییرات کارمایه پرداخته شده است. 
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شکل ٩‏ نمودار تغییرات میانگین خطای کارمایه تقریبی و میانگین خطای کارمایه دقیق صفحه نامحدود سوراخدار 
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لف) میانگین خطای دقیق ب) میانگین خطای تقریبی 


شکل ۱۰ توزیع میانگین حطای دقیق و تقریبی صفحه نامحدود سوراخدار 
در این نمونه رفتار براورد کننده خطا. در براورد روش ایزوژئومتریک بررسی‌می‌شود. 
حطای حاصل از تحلیل مسأله دارای یک نقطه تکین با بیشتر مسأله‌های با نقطه تکین؛ به علت بالا بودن 
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خطای آلودگی. مسأله‌های خحوبی برای بررسی کارایی 
برآورد کننده‌های خطا به شمار نمی‌روند. ولی هدف از 
این نمونه همان گونه که بحث شد. بررسی رفتار برآورد 
کننده خطا می‌باشد. نوک ترک به عنوان یک نقطه تکین. 
نسبت به دیگر نقطه‌ها است. خطای زیادی دارد. بنابراین 
یک برآورد کنندة حطای مناسب باید به نسبت بتواند این 
توزیم خطا را نمایش دهد. اگر بهبود شبکه با توجه به 
نتیجه‌های حاصل از برآورد کنندةءٌ خطا مورد نظر باشد. 
برآورد کنندهٌ حطایی که از آرایشی مناسب برخحوردار 
است. می تواند برنامه را در هر مرحله به بهبود محلی 
شبکه در نزدیکی نوک ترک راهنمایی کند. در شکل (۱۱) 
ویژگی‌های صفحه مربعی با تنش کششی ‏ مشاهده 
می‌شود. ضلع‌های این مربع به طول 28 و گسترش ترک 


به میزان 2 فرض شده است. 


ار تا ربق کر 
4 ۷۷۱۷۱۷۷ ۱۷ 


شکل ۱۱ صفحه مربعی ترکدار زیر اثر تتش کششی قائم 


بازشد گی ها در مود اول وت افش ان 
مختصات قطبی ی در شکل (۱۱) نتبجه‌های 


تحلیلی آن به ضورت زیر می‌باشد [20]. 


)۲۸( هت ات ِ 99 
2 2 2 93 

(۲۹( تهب ]هه ت 6 
2 2 2 مر ۷ 


سال بیست و دوم شماره دی ۱۳۹۰ 


۳۷ 


(۳۲۰) و 1 5 
2 .2 رل 9 
وراه تشه یت تشه ی تا رات 
(۳۱) به دست می‌آید. 
(۳۱ ۷۵ 2 ,1 

ات ارف ناشیا تعاها : 
الگوسازی می‌شود. عامل‌های به کار رفته در حل مسأله به 
صورت 10-< ظ : 1-1000 : 0.3 < ۷: 210 20 < طول 
و 

مسأله با روش ایزوژئومتریک برای تابع‌های شکل 
درجه یک. دو و سه تحلیل شده است. 

برای بررسی اثر افرايش درجه تابع‌های شکل در 
کارایی برآورد کننده خطاء شمار اجزا (بردار گرهی) در 
هر سه الگو ثابت فرض می‌شود بنابراین با توجه به رابطه 
۲1 با افزایش درجه تابع‌های شکل (۳) و نظر به 
اینکه شمار تکرار عضوهای صفر و یک بردار گرهی. 
1 می‌باشد. باید درهر جهت یکی به شمار نقطه‌های 
هادی (0) افروده شود. شکل (۱۲) بردار گرهی در هر 
ناحیه و اجزای تولید شده بر روی کل دامنه الگوسازی 
شده را با تابع‌های شکل درجه یک دو و سه نشان 
می‌دهد. 

و ها وتا هت وس 
شرایط مرزی و نقطه‌های هادی به کار رفته برای تحلیل 
صفحه ترکدار در کشش با تابع‌های شکل درجه یک دو 
ان کم تا اه( شا 
هادی در نزدیکی نوک ترک که خطای بیشتری نسبت به 
دیگر نقطه‌ها دارد» افزایش يافته اند. 

در شکل را ها وت الم ات که 
خا امه دق و رو را فا اه هن کزان 
در کشش با تابع‌های شکل درجه یک. دو و سه به 


نمایش درآمده ان 
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شکل ۲ (الف) بردار گرهی و (ب) اجزای تولید شده بر روی دامنه الگوسازی شده با تابع‌های شکل درجه یک دو و سه به روش ایزوژئومتریک 


نوک ترک 


ناحیه اول 
چپ پ چپ چپ پ پ پ چ چپ 
(الف) 
3 ۳9 
+ + چ ۷ ب ۲ 
(ب) (ج) 


شکل ۱۳ دامنه الگوسازی شده همراه با شرایط مرزی و نقطه‌های هادی به کار رفته برای تحلیل صفحه ترکدار در کشش با تابع‌های شکل (الف) 


درجه یک (ب) درجه دو (ج) درجه سه 
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1۹ 


میانگین خطای کارمایه دقیق 


ج‌( تابع‌های شکل درجه سه 


شکل ۱۶ مقایسه میانگین خطای کارمایه دقیق و تقریبی صفحه ترکدار در کشش با تابع‌ها شکل درجه یک دو و سه 


امین مگیم قطای کارسایه تترییی بزای هر 
سه گونه تابع‌های درجه یک دو و سه در شکل (۱۶) 
مناسب برای برآورد خطای صفحه ترکدار در کشش دارد. 
به گونه‌ای که در هر سه مورد میزان خطا در نوک ترک 
نسبت به دیگر نقطه‌ها به طور قابل ملاحظه‌ای بیشتر 
است. همچنین اگر مقایسه‌ای بین میزان بیشینه میانگین 
خحطای کنو دو شتسعل (۱۶) صسورت کت رگن مشاهده 
می‌شود که با افزایش درجه تابع‌های شکل خطای تحلیل 
ایزوژئومتریک کاهش می‌یابد. 


کز ارت مقاله حگونگی به کارگیری زوشی برای برآورد 
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خطای موجود در تحلیل مسأله‌ها به روش ایزوژئومتریک 
شرح داده شد. ویژگی بارز این روش هماهنگی آن با 
روش ایزوژئومتریک» که یک روش هندسی در تحلیل 
مسأله‌ها با معادلات دیفرانسیل جزئی حاکم مانند مسائل 
کشسانی است. به شمار می‌رود. با توجه به نتیجه‌های 
تحلیل برای شاخص اثر گذاری و همچنین برابری تقریبی 
میانگین خطای کارمایه دقیق و تقریبی برای دو مساله 
نمونه در این مقاله. می‌توان از ایين روش به عنوان 
راهکاری ساده و مهندسی برای برآورد خطا و بهبود 
میدان تنش به دست آمده از تحلیل مسأله‌ها به روش 


ایزوژئومتریک بهره گرفت. 
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